Electronic structure of normal and inverse spinel ferrites from first principles 
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We apply the self-interaction corrected local spin density approximation to study the electronic 
structure and magnetic properties of the spinel ferrites MnFe2 04, Fe3C>4, CoFe2 04, and NiFe2 04. 
We concentrate on establishing the nominal valence of the transition metal elements and the ground 
state structure, based on the study of various valence scenarios for both the inverse and normal 
spinel structures for all the systems. For both structures we find all the studied compounds to be 
insulating, but with smaller gaps in the normal spinel scenario. On the contrary, the calculated 
spin magnetic moments and the exchange splitting of the conduction bands are seen to increase 
dramatically when moving from the inverse spinel structure to the normal spinel kind. We find 
substantial orbital moments for NiFe2 04 and CoFe2C>4. 
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I. INTRODUCTION 

The field of spintronics is concerned with search for 
highly spin-polarized materials. One aim is to enhance 
tunnelling magnetoresistance (TMR) of magnetic tunnel 
junctions (MTJs) which are active members of magnetic 
random access memory (MRAM) elements. Also, the 
highly spin-polarized materials are of paramount impor- 
tance for increasing spin-polarization of currents injected 
into semiconductors, required for an optimal operation 
of spintronics devicesi There appears to be a number 
of ways to achieve high spin-polarization, most notably 
by employing fully spin-polarized ferromagnetic metals, 
namely half- metals (HM). 2 Another possibility is to ex- 
ploit features of the band structure of such tunnel bar- 
rier materials as e.g. MgO, and filtering electronic wave 
functions according to their symmetry in order to se- 
lect the most highly spin-polarized onesA 4 - The least ex- 
plored possibility is exploiting the spin-filtering effect, 
based on ferromagnetic or ferrimagnetic insulating barri- 
ers. It was introduced by Moodera et al£, using EuS tun- 
nel barriers. Spin filtering effect has been demonstrated 
in Gd/EuS/Al junctions^, which exhibit high magnetore- 
sistance, but show no great prospects for technological 
applications, on account of the low Curie temperature, 
T c , of EuS. 

Spinel ferrites 7 " have been studied for many years both 
regarding their magnetic behaviour and correlated na- 
ture in conjunction with their structural properties to in- 
crease their performance in high-frequency devices. Some 
of them can probably be used as spin filters. Spin- 
dependent gap should result in spin-dependent barrier 
for tunnelling of electrons through the insulator, giv- 
ing rise to spin filtering. Since the tunnelling proba- 
bility depends exponentially on the barrier height, the 
spin filtering efficiency can be very high. Candidates 
for spin filters include such spinel ferrites as NiFe204, 
CoFe204 and MnFe204». In particular, recently a spin 



filtering efficiency of up to 22% by the NiFe204 bar- 
rier has been reported by Liiders et alAifi In addition, 
Liiders et slM^ have demonstrated TMR of 120% in 
La2/3Sr 1 / 3 Mn03/SrTi03/NiFe204 junctions, which cor- 
responds to 45% spin polarization for the conductive 
NiFe204 film, which stays constant up to about 300 K. 

These spinel ferrites belong to the same family as mag- 
netite (Fe3 04) which has been most thoroughly stud- 
ied both for its HM character and the famous charge 
orde»iii4. Many theoretical studies have been dedi- 
cated to magnetite, and in particular its charge order 
(which will not be discussed in the present paper) us- 
ing various approximations to density functional the- 
ory (DFT) such as local spin density (LSD) approxima- 
tion or generalized gradient approximation (GGA), as 
well as going beyond, for example by invoking the Hub- 
bard U through the LDA+U (local density approxima- 
tion + U) approach, or using the self-interaction cor- 
rected (SIC)-LSD methodi 15 i 16 i 17 i 18 i 19 i 20 For the other 
spinel ferrites most theoretical studies have been done 
with LSD, GGA 8 ! 21 . 22 ! 23 or hybrid density functional a 4 . 
The former two approaches usually describe these ma- 
terials to be half-metallic and not insulating, if no dis- 
tortions are included. The reason is that the transition 
metal (TM) d electrons in oxides (as well as / elec- 
trons in rare earth compounds) are strongly correlated 
and cannot be adequately described within the stan- 
dard band theory framework with such approximations 
as LSD or GGA, placing them too high in energy around 
the Fermi level. The SIC-LSD method^ 5 ., on the other 
hand, provides better description of correlations than 
LSD, and was successfully applied to a variety of d— and 
/—electron materials^SLS. In this paper we apply the 
self-interaction corrected local spin density (SIC-LSD) 
approximation to study the electronic structure of spinel 
transition metal oxides MnFe204, Fe304, CoFe204, and 
NiFe2 04. We concentrate on the nominal valence of the 
TM elements and electronic and magnetic properties of 
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these systems in both normal and inverse spinel struc- 
tures. The reason being that in most cases these materi- 
als appear to exist as some mixture of those structures. 

The paper is organized as follows. An overview of the 
basic features of the SIC-LSD formalism is presented 
in the next section. Section III gives some computa- 
tional details, whilst the results of the application of the 
SIC-LSD method to the spinel ferrites TMFe 2 4 (with 
TM=Mn, Fe, Co and Ni), are presented and discussed in 
section IV. The paper is concluded in section V. 

II. THEORY 

The basis of the SIC-LSD formalism is a self- 
interaction free total energy functional, E SIC , obtained 
by subtracting from the LSD total energy functional, 
E LSD , a spurious self- interaction of each occupied elec- 
tron state tpc^, namely 

occ. 

E SIC = E LSD -J2€ W - (1) 
a 

Here a numbers the occupied states and the self- 
interaction correction for the state a is 

5% IC = U[n a ] + Ej£ D [n a ], (2) 

with J7[n a ] being the Hartree energy and E^ D [n a ] 
the LSD exchange-correlation energy for the correspond- 
ing charge density n a and spin density n a . It is the 
LSD approximation to the exact exchange-correlation 
energy functional which gives rise to the spurious self- 
interaction. The exact exchange-correlation energy E xc 
has the property that for any single electron spin density, 
n a , it cancels exactly the Hartree energy, thus 

U[n a ] + E xc [n a ] =Q. (3) 

In the LSD approximation this cancellation does not take 
place, and for well localized states the above sum can be 
substantially different than zero. For extended states in 
periodic solids the self-interaction vanishes. 

The SIC-LSD approach can be viewed as an extension 
of LSD in the sense that the self-interaction correction is 
only finite for spatially localized states, while for Bloch- 
like single-particle states E SIC is equal to E LSD . Thus, 
the LSD minimum is also a local minimum of E SIC . A 
question now arises, whether there exist other competi- 
tive minima, corresponding to a finite number of localized 
states, which could benefit from the self-interaction term 
without loosing too much of the energy associated with 
band formation. This is often the case for rather well lo- 
calized electrons like the 3d electrons in transition metal 
oxides or the 4/ electrons in rare earth compounds. It 
follows from minimization of Eq. that within the SIC- 
LSD approach such localized electrons move in a different 
potential than the delocalized valence electrons which re- 
spond to the effective LSD potential. Thus, by including 



an explicit energy contribution for an electron to localize, 
the ab initio SIC-LSD describes both localized and delo- 
calized electrons on an equal footing, leading to a greatly 
improved description of correlation effects over the LSD 
approximation, as well as, to determination of valence. 

In order to make the connection between valence and 
localization more explicit, it is useful to define the nom- 
inal valence^ as 

N va i — Z — N core — NsiC, 

where Z is the atomic number (26 for Fe), N core is the 
number of core (and semi-core) electrons (18 for Fe), and 
Nsic is the number of localized, i.e., self-interaction cor- 
rected, states (either five or six, respectively for Fe 3+ and 
Fe 2+ ). Thus, in this formulation the valence is equal to 
the integer number of electrons available for band for- 
mation. The localized electrons do not participate in 
bonding. To find the nominal valence we assume various 
atomic configurations, consisting of different numbers of 
localized states, and minimize the SIC-LSD energy func- 
tional of Eq. (JTJ with respect to the number of localized 
electrons. 

The SIC-LSD formalism is governed by the energetics 
due to the fact that for each orbital the SIC differenti- 
ates between the energy gain due to hybridization of the 
orbital with the valence bands and the energy gain upon 
localization of the orbital. Whichever wins determines 
if the orbital is part of the valence band or not, and in 
this manner also leads to the evaluation of the valence 
of elements involved. The SIC depends on the choice of 
orbitals and its value can differ substantially as a result 
of this. Therefore, one has to be guided by the energet- 
ics in defining the most optimally localized orbitals to 
determine the absolute energy minimum of the SIC-LSD 
energy functional. The advantage of the SIC-LSD for- 
malism is that for such systems as transition metal oxides 
or rare earth compounds the lowest energy solution will 
describe the situation where some single-electron states 
may not be of Bloch-like form. Specifically, in oxides, 
Mn-, Co-, Ni-, and Fe-3d states may be assumed to be 
localized, but not the O 2p states, because treating them 
as localized is energetically unfavourable. 

The SIC-LSD approach has been implemented^ 
within the linear muffin-tin-orbital (LMTO) atomic 
sphere approximation (ASA) band structure method 30 , 
in the tight-binding representational. In this method 
the polyhedral Wigner Seitz cell is approximated by 
slightly overlapping atom centered spheres, with a total 
volume equal to the actual crystal volume, while the 
electron wave functions are expanded in terms of the 
screened muffin-tin orbitals, and the minimization of 
E SIC becomes non-linear in the expansion coefficients. 
The so-called combined correction term 3 ^ has been 
implemented and consistently applied to improve on the 
ASA. 
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TABLE I: The lattice constants (a) and corresponding ASA 
radii (tasa) for the transition metal elements occupying tetra- 
hedral (tet) and octahedral [pet) sites, and for oxygen ions (in 
atomic units) for all the studied spinel ferrites. The radii of 
the empty spheres used in the calculations are not given, al- 
though typically four to five different types were used, all of 
the order of 1.65 - 1.95 atomic units, depending on the system, 
however, never exceeding the size of the oxygen spheres. 



System 


a 


tet 

?ASA 


„ oct 

?asa 


Ol 

?ASA 


02 

IASA 


MnFe 2 4 


16.08 


2.204 


2.588 


1.979 


1.979 


Fe 3 4 


15.87 


2.329 


2.752 


1.848 


1.848 


CoFe 2 4 


15.84 


2.434 


2.696 


1.928 


1.928 


NiFc 2 4 


15.78 


2.425 


2.686 


1.920 


1.920 



III. COMPUTATIONAL DETAILS 

The spinel ferrites of interest to the present study have 
a general chemical formula of the form AB2O4 and crys- 
tallize in the face-centred cubic structure. In the normal 
spinel structure A is a divalent element atom, occupying 
tetrahedral A sites, while B is a trivalent element, sitting 
on the octahedral B sites. When A is a trivalent element, 
and B consists of equal numbers of divalent and trivalent 
elements, distributed over crystalographically equivalent 
Bl and B2 octahedral sites, then the spinel structure is 
referred to as the inverse kind. In TM ferrites substantial 
off-stoichiometry and intersite disorder are often present 
in samples, but are not considered in this paper. The 
high temperature phase of magnetite is known to have 
the inverse spinel structure, where A atoms are Fe 3+ ions, 
and B sites are equally populated by Fe 2+ and Fe 3+ ions. 
Similarly, the NiFe2C>4 system has been established ex- 
perimentally to crystallize in the inverse spinel structure, 
with the A sites being Fe 3+ ions, while B sites equally 
populated by Ni 2+ and Fe 3+ ions. The MnFe2 04, on the 
other hand, is considered to be predominantly of the nor- 
mal spinel kind, as about 80% of A sites are populated 
by Mn 2+ ionsA^ 3 . The CoFe204 material is considered to 
be mostly an inverse spinel compound with about 80% 
of divalent Co ions occupying octahedral sites li 3 ^ 3 ^ As 
the experimental situation with respect to the observed 
structures and TM valences is not fully established, and 
the computer simulation of the exact physical conditions 
is very difficult, in this paper we study both extremes, 
namely the normal and inverse spinel structures for all 
the systems. In addition, we investigate a number of dif- 
ferent valence scenarios, defined in the following sections, 
to find the most energetically favourable solutions, be it 
only at zero temperature, for all the studied systems. 

Regarding the magnetic structure of the ferrites, we 
assume that of magnetite, with the spins of the TM atoms 
on the tetrahedral sublattice being antiparallel to those 
of the octahedral sublattice. Within a given sublattice 
the spins of all the TM atoms are arranged in parallel to 
one another. 

The calculations have been performed for the exper- 



imental lattice parameters, whose values, together with 
the corresponding ASA radii, are given in Table 0^. For 
the basis functions, we have used s-, p-, and d-muffin-tin 
orbitals on all the transition metal atoms as the so-called 
low waves and on the oxygen the s- and p-orbitals have 
been treated as low waves and the d-orbitals have been 
downfolded 3 ^. For a better space filling and to increase 
the number of basis functions, a set of empty spheres has 
also been included in the calculations. For the empty 
spheres only the s basis functions have been treated as 
low waves, while both p- and d-orbitals have been down- 
folded. All the calculations have been performed in the 
scalar-rclativistic mode, but for the calculated ground 
state configurations the spin-orbit coupling (SOC) was 
also included to calculate the orbital moment, in addi- 
tion to the spin moment. 



IV. RESULTS AND DISCUSSION 
A. MnFe 2 4 

As mentioned earlier, this compound is believed to be 
of predominantly normal spinel character. It is insulat- 
ing, with a small gap of 0.04-0.06 eV as determined by 
transport experiments^ Our calculations have adressed 
the important issues of this system by realizing both nor- 
mal ('N') and inverse (T) spinel arrangements of ions on 
the tetrahedral and octahedral sites. In addition, the all 
3+ scenario, where all the sites are occupied exclusively 
by the 3+ ions, has also been studied. Note that in the 
normal spinel environment the latter would mean that 
Mn 3+ (four d electrons are considered as localized) ions 
occupy the tetrahedral sites, while all the octahedral sites 
are exclusively populated by Fe 3+ ions ('N3+' scenario). 
For the inverse spinel environment the tetrahedral ions 
would be of Fe 3+ type, with the Bl octahedral sites oc- 
cuped by Mn 3+ ions and the B2 sites by Fe 3+ ions ('I3+' 
scenario) . In Table [HI we summarize the total energy 
differences for all the scenarios studied for MnFe2C>4, in 
comparison with all the other spinel ferrites. 

We find the normal spinel arrangement of ions to be the 
calculated ground state for MnFe2C>4 ( Table Hljl . in agree- 
ment with the experimental evidence for predominantly 
normal spinel character of this compound. The ground 
state solution is followed closely by the all 3+ scenario, 
realized in the normal spinel environment ('N3+'), which 
lies only 0.28 eV higher in energy. The inverse spinel so- 
lution is 0.58 eV higher, while the 'I3+' scenario is the 
most unfavourable state for MnFe2 04. 

To establish whether the degree of localization of d 
electrons of the TM ions residing on the tetrahedral sites 
has any bearing on the preference towards normal spinel 
structure in MnFe2 04, we have looked at the change in 
the localization energy when switching Mn 2+ ion between 
the tetrahedral and octahedral sites. We have found 
that the localization energy associated with Mn 2+ on the 
tetrahedral sites is 0.15 eV smaller than when Mn 2+ ions 
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TABLE II: Total energy differences (in eV per formula unit), 
calculated within SIC-LSD, between the ground state config- 
uration and other valence/structure scenarios for all studied 
spinel ferrites at the experimental lattice constant. The row 
marked by 'N' means normal spinel arrangement, where the 
tetrahedral sites are occupied by the divalent ions, namely 
Mn 2+ in MnFe 2 4 compound, Fe 2+ in Fe 3 4 , Co 2+ in 
CoFe 2 4 and finally Ni 2+ in NiFe 2 4 , while the octahedral 
sites are populated exclusively by Fe 3+ ions. Similarly, the 
row marked by T means that the Bl sites are occupied by 
Mn 2+ in MnFe 2 4 , Fe 2+ in Fe 3 4 , Co 2+ in CoFe 2 4 and 
Ni 2+ in NiFe 2 4 , with all the tetrahedral and B2 octahedral 
sites taken by the Fe 3+ ions. The notation 'I3+' means that 
tetrahedral sites and B2 octahedral sites in all the compounds 
are occupied by Fe 3+ ions, while the Bf octahedral sites are 
populated by Mn 3+ ions in MnFe 2 4 , Fe 3+ ions in Fe 3 4 , 
Co 3+ ions in CoFe 2 4 , and Ni 3+ ions in NiFe 2 4 . In the 
'N3+' scenario all the octahedral sites, in all the compounds 
studied, are occupied by the Fe 3+ ions, and the tetrahedral 
sites are taken by Mn 3+ ions in MnFe 2 4 , Fe 3+ ions in Fe 3 4 , 
Co 3+ ions in CoFe 2 4 , and Ni 3+ ions in NiFe 2 4 . Note, that 
in the Fe30 4 case the latter two scenarios, namely 'N3+' and 
'13+', are equivalent. 



MnFe 2 04, normal spinel 



Scenario MnFe 2 4 


Fe 3 4 


CoFe 2 4 


NiFe 2 4 


I 


0.58 


1.54 


0.20 


0.00 


13+ 


0.92 


0.00 


0.00 


0.52 


N 


0.00 


2.46 


1.09 


1.66 


N3+ 


0.28 


0.00 


0.46 


1.57 



occupy the octahedral sites. At the same time the local- 
ization energy associated with the Fe 3+ ions is smaller by 
0.19 eV for the tetrahedral sites, in comparison with the 
situation when on the octahedral sites. Hence, the local- 
ization energy alone favours the normal spinel structure 
only by 0.04 eV over the inverse spinel kind, which con- 
stitutes just a tiny fraction of the total energy difference 
of 0.58 eV (Table EI). Thus, the preference of MnFe 2 4 
for the normal spinel structure is mostly driven by other 
electronic degrees of freedom. 

Regarding the density of states (DOS) of MnFe2C>4 
in the calculated ground state normal spinel structure 
(the top panel of Fig. one can see that it is insu- 
lating, with a gap of about 0.075 eV which is slightly 
larger than the experimental value obtained from trans- 
port measurements. A larger gap of about 0.3 eV is seen 
for the inverse spinel structure (bottom panel), but the 
spin splitting of the conduction band is here considerably 
smaller than for the ground state normal spinel scenario 
(Table ITTfll . In variance to the normal and inverse spinel 
cases, the 'N3+' scenario (middle panel) is found to be 
half-metallic. 

For all the studied scenarios, the valence band is pre- 
dominantly of the oxygen type, with a very small ad- 
mixture of the TM character, and its polarization at the 
top of the band changes from positive to negative be- 
tween the inverse and normal spinel scenarios ( Table ITTT|l . 
What also changes substantially when moving from the 
normal to inverse spinel scenario is the conduction band 




E-E F (eV) 

MnFeiC>4, N3+ scenario 




E - E F (eV) 
MnFe-,0 4 , inverse spinel 




E - E F (eV) 

FIG. 1: (Color online) Spin decomposed total densities of 
states (in red), per formula unit, for MnFe 2 4 for the normal 
spinel structure ('N') (top), for the all 3+ scenario in the nor- 
mal spinel arrangement of atoms ('N3+') (middle), and for 
the inverse spinel scenario (T) (bottom). The oxygen con- 
tribution to the total density of states is also shown (green, 
dotted lines). The minority DOS is printed on the negative 
side of the y-axis, while the majority contributions are pre- 
sented on the positive side of this axis. 



splitting, which for the normal spinel scenario is about 
2.5 eV larger than that of the inverse spinel case. The 
reason being that the unoccupied TM (i-states are sub- 
stantially pushed up in energy when in the normal spinel 
environment. Considering that in reality MnFe2 04 is not 
a pure normal spinel compound, the exchange splitting 
of the conduction bands is most likely to be somewhere 
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TABLE III: Spin decomposed exchange splittings of the va- 
lence and conduction bands, as well as the energy gaps 
(in eV), for both inverse and normal spinel structures for 
MnFe2 04. Here VBM stands for the valence band maximum 
and CBM for the conduction band minimum, and f refers to 
spin-up- and J. to spin-down-component. 





T scenario 


'N' scenario 


VBM T - VBM 1 


0.69 


-0.91 


CBM T - CBM^ 


1.31 


3.85 


CBM T - VBM T 


1.64 


4.84 


CBM 1 - VBM 1 


1.02 


0.075 


Gap 


0.33 


0.075 



in between the values calculated for the normal and in- 
verse spinel scenarios. Of course, the larger the splitting, 
the more advantageous should it be for the spin filtering 
properties. 

To understand details of the densities of states shown 
in Fig. associated with the localized TM states, one 
should keep in mind that in the normal spinel scenario 
(top panel) which, as already mentioned, is the calculated 
ground state structure for this compound, the tetrahedral 
sublattice is populated by Mn 2+ ions and that in this case 
five minority Mn d electrons are described as localized 
states, seen just below the predominantly O 2p valence 
band, at about -6 eV. The unoccupied, majority Mn d 
electrons give rise to the states seen as peaks above the 
Fermi energy at about +6 eV. The octahedral sublattice 
is populated by Fe 3+ ions, with five majority d electrons 
localized, as the two sublattices are anti-parallel to one 
another. As a result the localized majority d Fe states 
are seen at about -12 eV below the Fermi energy. The 
unoccupied Fe d states, the minority ones, are seen just 
above the Fermi energy, over the range of up to about 
3.5 eV. 

In the 'N3+' scenario (middle panel), one minority Mn 
d electron gets delocalized to realize Mn 3+ ions on the 
tetrahedral sites. As a result the localized, minority Mn 
d peak has moved down in energy, lying just above -10 
eV, while the fifth, now delocalized, minority d electron 
appears at the Fermi energy, bringing down also the un- 
occupied, minority, Fe d states. The situation changes in 
the inverse spinel scenario (bottom panel), as now Mn 2+ 
ions reside on Bl sites of the octahedral sublattice, while 
the Fe 3+ ions populate the tetrahedral sublattice and B2 
sites of the octahedral sublattice. Thus on the tetrahe- 
dral sites we have five minority Fe d electrons which give 
rise to the peak at about -12 eV, while the localized Fe d 
electrons of the octahedral sublattice are seen as the peak 
on the majority side, at about -14 eV. The five localized 
Mn majority d states are again seen just below the va- 
lence band at about -7 eV. The unoccupied Mn minority 
d bands are seen at about 6 eV above the Fermi energy. 
The five unoccupied, majority, Fe d states, associated 
with the tetrahedral sublattice, are seen at about 2.5 eV 
above the Fermi energy. Finally, the unoccupied, mfnor- 



TABLE IV: Total spin magnetic moments (in /i_b per formula 
unit), calculated within SIC-LSD, for all the studied spinel 
ferrites and scenarios. 



Scenario MnFe2 04 


Fe 3 4 


CoFe 2 4 


NiFe 2 4 


I 


5.00 


4.00 


3.00 


2.00 


13+ 


4.10 


4.00 


3.00 


2.00 


N 


5.00 


6.00 


7.00 


8.00 


N3+ 


5.00 


4.00 


5.60 


6.80 



ity, Fe d states of the octahedral sublattice are seen as 
two separate peaks above the Fermi energy. Note, how- 
ever, that the SIC-LSD eigenvalues have no direct phys- 
ical interpretations as removal energies, and thus should 
not be directly compared with spectroscopies. To do so, 
one would need to take into account relaxation/screening 
effects that are not included in such an effective one- 
electron theory as SIC-LSD. One way to accomplish this 
is to employ the Ascf calculations, 38 ' 39,40 and another is 
the SIC-LSD based optimized effective potential (OEP) 
methodi 4n i 43 i 44 i 45 

The magnetic properties change when moving from 
the normal to inverse spinel scenario, as seen in Table 
IIVI where we compare the total spin magnetic moments 
for all the studied spinel ferrites. The total spin mag- 
netic moment for MnFe2 04 is 5 fj,B per formula unit, for 
both insulating and half-metallic solutions, while for the 
metallic '13+ ' scenario the spin moment is reduced to 
4.1 fis per formula unit. Note that unlike in the other 
ferrites there is no change in the total spin magnetic mo- 
ment between the normal and inverse spinel scenarios. 
The reason being that, as seen in Tabled there are only 
very small changes in the values of the spin moments 
of the transition metal ions, that are compensated by 
changes in the induced oxygen spin moments. 

TABLE V: Type-decomposed spin magnetic moments (in jis 
per formula unit), calculated within SIC-LSD, for MnFe 2 4 
for inverse and normal spinel structures. Here 'A' marks the 
tetrahedral-, while 'Bl' and 'B2' the octahedral-sites, and 
'01' and '02' stand for two different types of oxygens. 



Scenario FeA 3+ 


Mn sl 2 + 


Fe S2 3 + 


01 02 


I -4.09 


4.58 


4.11 


0.12 0.03 


Scenario Mn^ 2+ 


Fe B1 3 + 


Fe S2 3 + 


01 02 


N -4.49 


4.11 


4.11 


0.34 0.34 



Including spin-orbit coupling for the ground state nor- 
mal spinel scenario for MnFe2 04 compound, we find no 
considerable orbital moments either on Mn (-0.0005 /is) 
or on Fe ions (0.019 /is). As a result, the total orbital 
moment is of the order of 0.045 hb per formula unit, 
while at the same time the total spin moment is changed 
from 5.0 [is to 4.9995 [is per formula unit. Also, even 
with SOC included, we still observe a small energy gap 
of 0.064 eV, which incidently is in good agreement with 
the transport experiments^ 
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B. Fe 3 4 

Based upon its high magnetoresistive properties, mag- 
netite is of interest for technological applications, as e.g. 
computer memory, magnetic recording, etc. Magnetite is 
thought to be half-metallic, with the highest known T c of 
860 K. At about Ty = 122 K it undergoes a transition to 
an insulating state, associated with some kind of charge 
order, setting in on the octahedral sites, and a distortion 
of the crystal structure from the inverse spinel cubic to 
monoclinioiiiiii&. Verwey argued that below the transi- 
tion temperature, TV, the Fe 3+ and Fc 2+ cations order in 
the alternate (001) planes, and interpreted this transition 
as an electron localization-delocalization transition^. 

In the earlier paper, we have studied three different 
types of charge order on the octahedral sites, both in the 
high temperature (cubic) and low temperature (mono- 
clinic) phases 20 . In this paper, for the sake of comparison 
with other spinel TM oxides, we concentrate exclusively 
on the high temperature cubic phase and the scenarios 
enumerated in Table llVl 

As seen from Table |n] and Fig. [2 (middle panel), 
the ground state of magnetite in the cubic phase is 
half-metallic, with all Fe ions in 3+ configuration (five 
d electrons localized). This ground state scenario 
(T3+'='N3-|-') is the result of a derealization of the sixth 
d electron of the original Fe 2+ ions, that together with 
Fe 3+ ions randomly populate the octahedral Bl and B2 
sites in the high-temperature cubic phase. In the ground 
state scenario, this sixth electron is seen to give rise to 
the peak at the Fermi energy in the minority channel, 
together with the other 10 unoccupied minority d states 
associated with the octahedral sites. The five localized 
tetrahedral Fe^ minority d states appear around -13 eV, 
while the unoccupied majority Fe,4 d states are seen just 
above the Fermi energy. All the localized majority d 
states of the octahedral sites are at about -12.5 eV. The 
valence band is of predominantly O p character. 

The inverse spinel solution (Table |n] and the bottom 
panel of Fig. [2J, corresponding to the assumed Verwey 
charge order, lies about 1.5 eV above the ground state. 
In this scenario, the above mentioned sixth electron is 
localized on the Bl-octahedral sites, and appears as a 
small hump at the bottom of the minority valence band. 
The remaining five d states of the Fesi sites are seen 
at about -8 eV below the Fermi energy in the majority 
bands. The five localized majority Fes2 d states lie at 
about -12.5 eV. The minority Fe^ d states (seen below 
-10 eV) are localized, and the majority Fe^ d states are 
unoccupied, occuring at about 3 eV above the Fermi en- 
ergy. As a result, since the four remaining minority Fe^i 
d states and five minority FeB2 states are also unoccu- 
pied, our calculations for this scenario give an insulating 
state with a gap of ~0.7 eV. Understandably, as the lat- 
ter has been calculated for the high-temperature inverse 
spinel structure, its value is much larger than the exper- 
imental value of 0.14 eV—, measured for the true, low 
temperature monoclinic phased 



Fe^Oa, normal spinel 
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FIG. 2: (Color online) Spin decomposed total densities of 
states (in red), per formula unit, for FesCU for the normal 
spinel structure ('N') (top), for the all 3+ scenario in the in- 
verse spinel arrangement of atoms ('13+') (middle), and for 
the inverse spinel scenario (T) (bottom). The oxygen con- 
tribution to the total density of states is also shown (green 
dotted). The minority DOS is printed on the negative side of 
the y-axis, while the majority contribution is shown on the 
positive side of this axis. 



Our calculations for the normal spinel structure give 
the most energetically unfavourable solution for mag- 
netite, lying about 2.5 eV above the T3+' ground state 
scenario. In this normal spinel case (top panel of Fig. |5J), 
we obtain an insulating solution with an energy gap of 
0.08 eV. Here, the tetrahedral sites are occupied by Fe 2+ 
ions, while the octahedral sites are populated with Fe 3+ 
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TABLE VI: Spin decomposed exchange splittings of the va- 
lence and conduction bands, as well as the energy gaps (in 
eV), for both inverse and normal spinel structures for FesCU. 
Here VBM stands for the valence band maximum and CBM 
for the conduction band minimum, and j refers to spin-up- 
and I to spin-down-component. 





T scenario 


'N' scenario 


VBM T - VBM 1 


0.06 


-0.35 


CBM T - CBM^ 


1.61 


3.65 


CBM T - VBM T 


2.33 


4.08 


CBM 1 - VBM 1 


0.78 


0.08 


Gap 


0.72 


0.08 



TABLE VII: Total and type-decomposed spin magnetic mo- 
ments (in Bohr magnetons per formula unit) for magnetite as 
calculated within SIC-LSD for three different scenarios. Only 
different Fe-types are listed in the table. 



Scenario M to tai 


M Fe 2+ 


M Fe 3+ 


M Fe 2 + 


M Fe 3+ 

m 




I 


4.00 




-4.00 


3.57 




4.08 


13+ 


4.00 




-4.02 




3.90 


3.90 


N 


6.00 


-3.46 






4.09 


4.09 



ions. As a result, five localized minority Fe^ d states 
lie at about -8 eV, while the sixth localized, majority, 
Fe,4 d state sits right at the bottom of the majority va- 
lence band. All the remaining unoccupied majority Fe^ d 
states are just about seen at 5 eV above the Fermi energy, 
giving rise to a large exchange splitting of the conduction 
band. The localized, majority, Fesi and Fes2 d states 
are seen at -12.5 eV, while their unoccupied, minority, 
states lie just above the Fermi energy, over the range of 
about 3-4 eV. So, again like in MnFe2C>4, we see a large 
change in the exchange splitting of the conduction band 
when moving from the inverse to normal spinel structure. 
For the inverse spinel scenario (T) it is of the order of 
1.6 eV, while for the normal spinel arrangement it in- 
creases to 3.65 eV (Table EB- It is interesting to note, 
that in the normal spinel structure the sixth localized d 
electron of the Fe 2+ ion occupies one of the e g orbitals, 
while in the inverse spinel scenario it populates one of 
the t2g states. 

The total spin magnetic moment per formula unit is 4 
[iB (Tabic ITV|) for all the scenarios studied, with the ex- 
ception of the normal spinel scenario, where we see a 50% 
increase to 6 [is- As all the scenarios give rise to either 
insulating or half-metallic states, the spin magnetic mo- 
ments are naturally integer numbers. Table IVTT1 explains 
how the 50% increase in the total spin magnetic moment 
comes about when switching from the inverse to normal 
spinel structure. In the inverse spinel case the spin mo- 
ment of the tetrahedral Fe ions gets just about cancelled 
by the spin moment of the B2-octahedral sites, so that 
the total spin moment is mostly due to the Fe 2+ ions on 
the Bl-sites. In the normal spinel, on the other hand, 
the spin moment of the tetrahedral Fe 2+ ions is smaller, 



and oppositely alligned with the spin moments of the 
Fe 3+ ions that occupy all the octahedral sites. Bearing 
in mind that there are twice as many octahedral sites 
as the tetrahedral ones, the substantial increase is easy 
to account for, especially that the induced spin moments 
on the oxygen sites do not differ much between the two 
scenarios. 

Including spin-orbit coupling for the ground state '13+ ' 
scenario leads to a very small total orbital moment of 
about 0.05 \ib per formula unit, while the total spin mo- 
ment is very slightly reduced from 4 \±b to 3.9998 \xb per 
formula unit. The orbital moments due to the individual 
Fe-ions are similarly very small, with the tetrahedral Fe 
being -0.015 \xb and the octahedral Fe of 0.035 \ib- 



C. CoFe 2 4 

This compound is believed to be mostly of the inverse 
spinel kind42i£2iiii, with divalent Co ions occupying pre- 
dominantly the octahedral sites. However, similarly to 
magnetite in the high temperature cubic phase, our cal- 
culations find the ground state of CoFe204 to be half- 
metallic and of the '13+ ' type (Table and Fig. 
middle panel). As seen in Table ITT1 the inverse spinel 
scenario (T) is not far, lying only 0.2 eV higher in en- 
ergy. What this seems to imply is that this compound 
prefers the inverse arrangement of atoms, independently 
of the actual valence of the Co ions. The normal spinel 
solution is about 1.1 eV away. Although the ground state 
we find is half-metallic, the inverse spinel scenario ('I'), 
with Co 2+ ions occupying the Bl octahedral positions, 
describes CoFe204 as an insulator, with a gap of 0.8 eV, 
which is reduced to 0.21 eV in the normal spinel case 
(Table IvTlTl and Fig. |3J . 

To understand in detail the densities of states in Fig. 
[3] note that Co has only one minority electron more than 
Fe. So, the Co 2+ ion has two localized minority d elec- 
trons, in addition to the five majority ones. In the ground 
state 'I3+' scenario (middle panel in Fig. |3J), one of these 
two minority electrons gets delocalized, contributing to 
the states at the Fermi energy, while the other, local- 
ized, one is seen as a sharp peak just above -10 eV. All 
the remaining main features of DOS for this scenario are 
exactly like in the case of 'I3+' scenario in magnetite. 

For the inverse spinel structure (bottom panel of Fig. 
|3J|, the situation is again very much like in magnetite 
(bottom panel of Fig. with the exception that now 
we have a small double hump, at about -7 eV, slightly de- 
tached from the predominantly O 2p valence band, while 
in magnetite it was still attached to the valence band and 
represented just a single minority d electron. 

In the normal spinel scenario the Co 2+ ions now reside 
on the tetrahedral sites, with their five localized minority 
d states seen as a rather sharp peak at about -8 eV, in the 
top panel of Fig. [3J The remaining two localized majority 
d states are sitting just below the majority valence band. 
All the other features are like in magnetite. 
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CoFe-,0.}. normal spinel 




CoFe-,0 4 , inverse spinel 





E-E F (eV) 

FIG. 3: (Color online) Spin decomposed total densities of 
states (in red), per formula unit, for CoFe204 for the normal 
spinel structure ('N') (top), for the all 3+ scenario in the in- 
verse spinel arrangement of atoms ('13+') (middle), and for 
the inverse spinel scenario (T) (bottom). The oxygen con- 
tribution to the total density of states is also shown (green 
dotted lines). The minority DOS is printed on the negative 
side of the y-axis, while the majority contribution is shown 
on the positive side of this axis. 



TABLE VIII: Spin decomposed exchange splittings of the 
valence and conduction bands, as well as the energy gaps 
(in eV), for both inverse and normal spinel structures for 
CoFe2 04. Here VBM stands for the valence band maximum 
and CBM for the conduction band minimum, and f refers to 
spin- up- and 1 to spin-down-component. 





T scenario 


'N' scenario 


VBM T - VBM 1 


0.22 


-0.24 


CBM T - CBM 1 


1.28 


4.07 


CBM T - VBM T 


2.08 


4.52 


CBM 1 - VBM 1 


1.02 


0.21 


Gap 


0.80 


0.21 



TABLE IX: Type-decomposed spin magnetic moments (in \xb 
per formula unit), calculated within SIC-LSD, for CoFe2 04 
for inverse and normal spinel scenarios. Here 'A' marks the 
tetrahedral-, while 'Bl' and 'B2' the octahedral-sites, and 
'01' and '02' stand for two different types of oxygens. 



Scenario FeA 3+ 


Co sl 2+ 


Fe S2 3+ 


Ol 


02 


I -4.11 


2.58 


4.11 


0.13 


0.07 


Scenario Coa 2+ 


Fe sl 3+ 


Fe S2 3+ 


Ol 


02 


N -2.58 


4.13 


4.13 


0.32 


0.32 



normal spinel structure, the e g minority states are pop- 
ulated before the t2 9 ones, which is in variance to the 
inverse spinel structure. In the latter case, the Co 2+ ions 
reside on the Bl octahedral sites, and the tz g states are 
lying lower in energy than the e g states. 



As seen in Table ITVl the total spin magnetic moment 
for both the T and '13+ ' scenarios is 3 /is per formula 
unit. It is reduced from 4 fis in magnetite due to the 
smaller value of the spin moment of the divalent Co-ion 
of 2.58 hb (Table ITxTl . in comparison with the spin mo- 
ment of the divalent Fe ions of 3.57 /ib fTabl dViTfr . What 
is however more dramatic is the change of the spin mo- 
ment when moving from the inverse to the normal spinel 
arrangement of ions. Table IIVI shows that in the nor- 
mal spinel scenario the total spin moment is 7 \xb per 
formula unit, which again is due to the fact that the 
octahedral sites are populated exclusively by Fe 3+ ions, 
whose spin moments are arranged in parallel to one an- 
other and whose magnitudes are considerably larger than 
the moment of Co 2+ ions on the tetrahedral sites. 



Similarly to other ferrites, one sees substantial change 
in the exchange splitting of the conduction band between 
the inverse and normal spinel scenarios, from 1.28 eV to 
4.07 eV. Also, the negative spin polarization of the va- 
lence band is seen in the normal spinel, while a positive 
one in the inverse spinel structure. Like in magnetite, 
for the divalent Co 2+ ions on the tetrahedral sites in the 



With respect to spin-orbit coupling we find the total 
orbital moment of the ground state '13+ ' scenario to be 
quite substantial of the order of 0.58 /is per formula unit, 
and associated mostly with the Co 3+ ion. As the total 
spin moment is reduced from 3 [Ib to 2.9997 /is per for- 
mula unit, the ratio of the total orbital to spin moment 
is 0.19. The ratio of the orbital to spin moment for the 
Co 3+ ion itself is 0.21. 
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D. NiFe 2 4 

NiFe2 O4 is a ferromagnetic insulator that is of possible 
interest as a spin filter in MTJs 9,10 . This compound has 
the Curie temperature of 850 K, and hence has a great 
potential for technological applications. 

In agreement with experiments, we find the ground 
state of NiFe204 to be insulating and of the inverse spinel 
kind (Table |IV| and the bottom panel of Fig. EJ}. Thc 
calculated energy gap is 0.98 eV ( Table IXjl . It gets re- 
duced to 0.26 eV in the normal spinel scenario, which 
is the most energetically unfavourable solution for this 
compound. In thc ground state T scenario, the octahe- 
dral Bl-sites are occupied by Ni 2+ and B2 sites by Fe 3+ 
ions, while the tetrahedral sites are populated exclusively 
by Fe 3+ ions. Replacing the Co 2+ ions on the octahedral 
sites in CoFe2C>4 by Ni 2+ leads to the reduction of the 
total spin magnetic moment of 3 \xb in CoFe204 to 2 /is 
in NiFe2C>4 fTable llV)l . since the spin magnetic moment 
of the Ni 2+ ion is 1.57 [is (Table IXlJ) . as compared to 
2.58 hb spin moment of the Co 2+ ion ( Table llXjl . Thc 
oxygen spin moments in both materials are comparable, 
and aligned in parallel to the cation spin moments on the 
octahedral sites. The width of the predominantly oxygen 
2p valence band in NiFe2 04 is comparable to the one of 
CoFe 2 04, but is reduced with respect to the valence band 
of magnetite. The reason being that the sixth localized 
d electron of the Fe 2+ ion (in Fig. (middle panel) seen 
at the bottom of the valence band between -5.0 and -6.0 
eV) is strongly hybridized with the oxygen p band. The 
situation is different in NiFe2 04, where the three local- 
ized minority t2 9 electrons, seen at about -8.0 eV (Fig. 
bottom panel), are well separated from the bottom of the 
valence band. Also, the exchange splitting of the conduc- 
tion band, of importance to spin filtering, is about 20% 
smaller in NiFe204 than in the Verwey phase of Fe304. 



TABLE X: Spin decomposed exchange splittings of the va- 
lence and conduction bands, as well as the energy gaps 
(in eV), for both inverse and normal spinel structures for 
NiFe204. Here VBM stands for the valence band maximum 
and CBM for the conduction band minimum, and f refers to 
spin-up- and J. to spin-down-component. 



NiFeoO^ normal spinel 





T scenario 


'N' scenario 


VBM T - VBM 1 


0.10 


-0.12 


CBM T - CBM 1 


1.21 


2.93 


CBM T - VBM T 


2.19 


3.31 


CBM 1 - VBM 1 


1.08 


0.26 


Gap 


0.98 


0.26 




NiFei0 4 , inverse spinel 




E-E F (eV) 

FIG. 4: (Color online) Spin decomposed total densities of 
states (in red), per formula unit, for NiFe20 4 in the nor- 
mal spinel structure ('N') (top), NiFe204 in the inverse spinel 
structure for all 3+ ('13+') scenario (middle), and NiFe204 in 
the inverse spinel structure (T) (bottom). The oxygen con- 
tribution to the total density of states is also shown (green 
dotted lines). The minority DOS is printed on the negative 
side of the y-axis, while the majority contribution is shown 
on the positive side of this axis. 



To understand details of the DOS of the 'I3+' sce- 
nario of NiFe2 04, it is helpful to follow the discussion 
of the same scenario for CoFe204, keeping in mind that 
a Ni 3+ ion has two minority d electrons in addition to 
the five majority ones, localized by SIC. Obviously, as 
already mentioned when discussing MnFe204, the po- 
sitions of the localized d peaks calculated in the SIC- 



LSD should not be directly compared with photoemission 
experiments^ 

The changes in the electronic structure of NiFe204 
when moving from the ground state inverse spinel struc- 
ture to the normal spinel scenario arc immediately obvi- 
ous from comparing the bottom and top panels of Fig. 01 
In the normal spinel case the Ni 2+ ions occupy the tetra- 
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TABLE XI: Type-decomposed spin magnetic moments (in \xb 
per formula unit), calculated within SIC-LSD, for NiFe2C>4 
for inverse and normal spinel scenarios. Here 'A' marks the 
tetrahedral-, while 'Bl' and 'B2' the octahedral-sites, and 
'01' and '02' stand for two different types of oxygens. 



Scenario Feyt 


Ni B1 2 + 


Fe S2 3+ 


Ol 


02 


I -4.11 


1.57 


4.11 


0.14 


0.07 


Scenario NLt 2+ 


Fe B1 3 + 


Fe S2 3+ 


Ol 


02 


N -1.65 


4.13 


4.13 


0.35 


0.35 



hedral sites, while the octahedral sites are solely taken 
by the Fe 3+ ions. As a result, the total spin magnetic 
moment is increased from 2 fj,s per formula unit to 8 /xs 
per formula unit, as seen in Table ITT| in agreement with 
experimental findings of Refs. Illlll2l However, as seen in 
the top panel of Fig. the density of states is still just in- 
sulating in both spin channels. Also, unlike in the case of 
the inverse spinel structure, the valence band is strongly 
spin-polarized, and the polarization is negative. The oxy- 
gen spin magnetic moment is 0.35 /is, aligned in parallel 
to the Fe spin moment (see Table lX*T|) . and three times the 
value it has in the inverse spinel structure. Also, the Ni 
spin moment is slightly increased in magnitude to -1.65 
/.t b ■ Moreover, the exchange splitting of the conduction 
band is more than twice increased in the normal spinel, 
in comparison with the inverse spinel structure, from 1.21 
eV to 2.93 eV. As in the case of the other spinel ferrites 
in realizing Ni 2+ ions, the e g states are populated first, i. 
e., are lying lower in energy than the t2 9 states, which is 
opposite to the inverse spinel structure. However, ener- 
getically, the normal spinel structure for NiFe2 04 is very 
unfavourable with respect to the inverse spinel structure 
(Table HU. 

Including the spin-orbit coupling for the ground state 
inverse spinel scenario gives rise to the total orbital 
moment of 0.67 (ib per formula unit, and is mostly 
due to Ni ions, with some minor contributions from Fe 
atoms. This calculated value is over two times larger 
than calculated from LSD in Ref. 22. Also, in the 
earlier SIC-LSD calculations for TM oxides Svane and 
Gunnarsson» obtained the orbital moment for NiO 
of 0.27 fiB, which is substantially smaller than in the 
present calculations. Since the total spin moment is 
slightly reduced to 1.9997 /ib per formula unit, when 
SOC is taken into account, hence the ratio of the total 
orbital and spin moments is calculated to be about 0.34. 
The latter is in good agreement with the experimental 
estimates of 0.27±0.07 5 - and 0.34 (within error bars of 
up to ±0.11) for Ni in NiFe 2 4 and NiO^S Note that 



the orbital moment due to the Ni 2+ ion alone is about 
0.7 /is, while its spin moment is 1.58 \jlb-, both per 
formula unit, giving rise to the orbital to spin moment 
ratio of 0.44 for this ion. The total orbital moment 
is mostly due to the Ni 2+ ion, as the contributions of 
other TM ions are smaller by an order of magnitude or so. 

V. CONCLUSIONS 

We have shown that, owing to a better treatment of 
correlations, SIC-LSD can provide useful insights to the 
nature of a number of spinel ferromagnetic insulators. 
We have been able to address the issues of the normal 
versus inverse spinel arrangements in these systems, their 
electronic and magnetic properties and the valence of the 
transition metal atoms. We find all the studied ferrites to 
be insulating for both the inverse and normal spinel sce- 
narios, with the calculated energy gaps being smaller in 
the normal spinel environment, however showing an in- 
creasing trend when moving from MnFe2 04 to NiFe2 04. 
We have observed dramatic increase of the calculated 
spin magnetic moments, as well as the exchange splitting 
of the conduction bands, when moving from the inverse 
to normal spinel scenarios, some of which have been ob- 
served in experiments^ 

The total energy considerations seem to favour the in- 
verse spinel arrangement of TM ions as the ground state 
configurations for all the studied ferrites, with the possi- 
ble exception of MnFe2 04 where the normal spinel envi- 
ronment may be preferred, with Mn 2+ ions on the tetra- 
hedral sites and the Fe 3+ ions on the octahedral sublat- 
tice. Also, based on the total energy arguments, we find 
a partial derealization of the minority spin states to be 
favourable in Fe3 04 and CoFe204, leading to the half- 
metallic ground states with all TM ions in the trivalent 
configuration and in the inverse spinel arrangement. The 
fully inverse spinel scenario, with the Fe 3+ ions on the 
tetrahedral sites and the octahedral sites occupied both 
by Ni 2+ ions and Fe 3+ ions, is found to be the ground 
state only in NiFe204. Finally, these findings constitute 
a good starting point for further studies, incorporating 
alloying of the normal and inverse spinel structures. 
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